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In the context of non-relativistic quantum mechanics, we obtain several upper and lower limits on 
the mean square radius applicable to systems composed by two-body bound by a central potential. 
A lower limit on the mean square radius is used to obtain a simple criteria for the occurrence of 
S-wave quantum halo sates. 
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I. INTRODUCTION 

(N ; 

The root mean square (rms) radius is used in many fields of physics to characterize the size of quantum systems; 
, this observable is thus of special interest. In the case of two-body systems, the rms radius and the energy of one 
' eigenstate determine the depth and the range of the central potential which binds the particles provided that the 
^vq _ number of bound states supported by the potential is also known. A simple example is the naive description of the 
deuteron by a square- well potential (see for example pj). The depth Vo and the range R can be adjusted to reproduce 
the rms radius and the binding energy of the deuteron. This adjustment is not unique excepted if the number of 
bound states in the potential is fixed to one. Obviously, these quantities are not sufficient to infer the shape of the 
potential. Indeed, the same simple description of the deuteron could be achieved with an exponential potential for 
example. However, if information on the shape of the potential is obtained by other means then constraints on the 
I ' rms radius, on the energy and on the number of bound states yield strong restrictions on the potential. Consequently, 
j— i \ upper and lower limits on these quantities are interesting tools to obtain easily constraints on the interaction. 
There exists a fairly large number of upper and lower limits on the energy of eigenstates in the literature 
0000 0] as well as on the number of bound states supported by central potentials 0, 0, 0, 0, Q, Il5l TigTfl^ 
Similar results concerning the rms radius are scarcer. A first general inequality gives a lower limit on the rms radius 
O of a £-wave state in terms of the average kinetic energy (h = 2m = 1) jlsL p. 73] 

S : /jx > ( 2 ^ + 3 ) 2 m 

{r ) > 4(T) ■ (1) 

^ ' This relation is however not very useful except in the case of power-law potentials, V(r) — sgn(p)gr p (p > —2), for 
J> , which the virial theorem gives a simple relation between the energy, E, of the £-wave state and the average kinetic 
' energy: (T) = pE/ (p + 2). In this case, the relation |QJ simply reads 

M\ (r 2 )>(2£ + 3) 2 ^. (2) 

In the case of a vanishing angular momentum, t = 0, another restriction is given by the Bertlmann-Martin inequality 

ID 

This simple relation is actually also applicable to systems composed by N > 2 identical particles provided that no 
symmetry is required for the wave function |20| . However, this formula yields restrictions only on the size of the 
ground state of the system and the energy of two levels needs to be known. 

In sec. [n] we propose several rigorous upper and lower limits on the rms radius applicable to systems composed 
by two-body bound by a central potential. Some of these limits are applicable to arbitrary eigenstates of the systems 
(any value of the radial quantum number, n, and of the angular momentum £). These limits involve the energy of the 
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eigenstate considered as well as the potential itself. Upper and lower limits on the energy of the eigenstate can then 
be used to obtain upper and lower limits on the rms radius as a function of the potential only and to infer constraints 
on the values of its parameters. In sec. IIIII as a simple application of the results obtained in this work, we use a 
lower bound on the rms radius to study weakly bound systems (quantum halo states) and to find a criteria for the 
occurrence of such states. Some tests of the other limits are reported in sec. IIVI At last, we present some conclusions 
in sec. 

II. UPPER AND LOWER LIMITS ON THE RMS RADIUS 

To obtain various upper and lower limits on the rms radius, we consider the Schrodinger equation with a central 
potential 



u"{r) = 



«('•), (4) 



where u(r) — rR(r) and R(r) is the radial wave function; the angular part is obviously given by the spherical harmonics 
Yi m {9, if). For simplicity of the notations, we do not write the indices n and £ on the function u(r) and on the energy 
E, except if this is necessary for the clarity of the formula. The negative part of potential is supposed to be less 
singular than r~ 2 at the origin and the potential is supposed to be piecewise continuous elsewhere. 

To obtain the relation from which each limit will be derived in this work, we multiply the relation Q by r 2 u(r) 
and we integrate to obtain 

/>oo />oo 

/ drr 2 u(r)u"(r) = drr 2 V(r) u(r) 2 + £(£ + 1) - E(r 2 ). (5) 
Jo Jo 

Integration by parts of the normalization condition of the wave function u(r), J °° dru(r) 2 — 1, leads to 

2 / drru(r)u'{r) = -1. (6) 



In all cases throughout this paper, it is easy to verify that the boundary terms of all integration by parts are vanishing 
if the energy E and the rms radius (r 2 ) have a sense. Integration by parts of the left hand side of the relation (JSJ) 
together with 10 leads to 



-E(r 2 ) +£(£+!) = 1- drr 2 V(r)u{ry - drr 2 u(ry. (7) 

Jo Jo 

With the help of the equality Q, it is possible to obtain various upper and lower limits on the rms radius. We 
gives these limits in the next sections. 



A. Simple upper limit 

The first upper limit we present is directly obtained from the relation Q. Since J °° drr 2 u'(r) 2 > 0, we have 

/>oo 

-E(r 2 ) < 1- / drr 2 V(r)u(r) 2 - £(£ + 1) 
Jo 

< 1+ sup [-r 2 V{r)] -£(£+ 1). (8) 

0<r<oo 

We still need to show that the right hand side of the last inequality of © is finite, to have non trivial results, and 
positive for non vanishing value of the angular momentum. This last inequality is finite only if the negative part of the 
potential decreases as r~ 2 at infinity or faster. It is well known that the class of potentials characterized by a negative 
part which decrease faster than r~ 2 at infinity support only a finite number of bound states |2lj |. Consequently, there 
exist a maximal value, L, of the angular momentum, £, above which no bound state exists. If L + is defined by the 
relation 

sup [~r 2 V(r)] =(l+ + 1) , (9) 

0<r<oo \ z / 
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then it is also well known that L + is an upper limit on L, L < L+ [H|2l|. The upper limit JSJ takes then the form 

~E(r 2 ) < j + {L+f - I 2 + L+ - £, (10) 

which proves that its right hand side is indeed positive. 

Potentials characterized by negative parts which decrease exactly as cr~ 2 at infinity are at the borderline between 
potentials which possess a finite number of bound states and potentials for which this number is infinite. It has been 
proved that these potentials support only a finite number of bound states if c is small enough 22] . Thus in this case, 
there exists a maximal value, L, of the angular momentum. The upper limit L + is expected to be valid and the 
relation (|lf)|> proves that the right hand side of ijBJ is positive. However, to be completely rigorous, one should study 
this class of potentials in detail and prove that indeed that the upper limit L + is correct. 

The simple upper limit ijHJ shows that the product of the absolute value of the energy by the mean square radius 
cannot increase faster than linearly with the depth of the potential. 

This upper limit will be significantly improved in the next section. We have presented this result because it is very 
simple and the discussion about the positivity of the right hand side of the inequality JSJ will be useful later. 



B. Main results 



The result JSJl obtained in the previous section can be improved provided we take into account the contribution of 
the last term of Eq. (jTJ). Such a contribution is easily obtained. An integration by parts leads to the equality 

) 2 /*°° 

drr 2 u'(r) 2 = -- / dr rV \r)u" (r). (11) 
3 Jo 

The second derivative of the wave function which appears in (|11|) can be replaced using the Schrodinger equation (0J . 
A second integration by parts yield the desired result 

dr r 2 u'(r) 2 = ~ J™ dr (r 3 V(r))' u(r) 2 + l M±R _ E (r 2 ). (12) 

The relation Q together with the identity 11 2|) leads to 

4 If 00 
-2E(r 2 ) + -£(£+1) = 1 + - drW(r)u{r) 2 , (13) 
3 3 J 

where 

W(r) = -(6V{r) + rV'(r))r 2 . (14) 
The upper limit on the rms radius simply reads 

-E(r 2 ) < 1 + i sup [W{r)\ - 2 -l{£ + 1). (15) 

This upper limit l|15fl is non trivial if the negative part of the potential decreases as r~ 2 at infinity or faster and if the 
positive part of the potential is less repulsive than r~ 6 at the origin. The positivity of the right hand side of the upper 
limit l)15|l is simple to prove. In sec. Ill Al we have obtained that the quantity sup [— r 2 F(r)] — 1(1 + 1) (see Eq. 
is positive. In contrast, we need now to show that the quantity sup [W^r - )] /4 — £(£ + 1) is positive. For this purpose 
it is sufficient to prove that the inequality sup [— r 2 y(r)] < sup [^(r)] /4 is verified. Suppose that the supremum of 
the function — r 2 V(r) is reached for r — f. At this point, we have 

-\r"V'{f) = \f 2 V{f). (16) 

This infers that value of the supremum of —r 2 V(r) is equal to the value of W(r)/A at this point: —f 2 V(f) = W(f)/A. 
Thus the two functions have a crossing at r = f implying that the supremum of W(r)/4 cannot be smaller than the 
supremum of — r 2 V(r). 

Similarly, the lower limit obtained from (|13JI reads 

-E(r 2 ) > \ + \ inf [W(r)\ - + 1). (17) 
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The inequality (|17|l is non trivial if the positive part of the potential decrease as r~ 2 at infinity or faster and is less 
singular than r~ 2 or more singular than r~ 6 at the origin. For potentials less singular than r~ 2 at the origin, we 
have W(0) = yielding inf[W(r)] < 0. Consequently, the lower limit is only non trivial for S-wave states in this case. 
Moreover, if the infimum of W(r) is negative, the lower limit becomes trivial for a depth of the potential large enough. 

A case where the lower limit (II Til leads certainly to non trivial results concerns S-wave states and potentials such 
as W(r) is non negative. In this case we obtain the simple lower limit 

-2E e=0 (r 2 ) e=0 > I- (18) 

This relation is used in sec. IHII to obtain a criteria for the occurrence of S-wave quantum halo states. We show that 
this criteria is still efficient even if the function W(r) associated to the potential is somewhere slightly negative. Note 
also that the lower limit (|18|l has been obtained previously but for a more restrictive class of potentials composed of 
potentials of finite range, V(r > R) = 0, which are less singular than r~ 2 at the origin [23l |. 

It is interesting to note that, contrary to Eq. JSJl, the upper and lower limits l(T5|) and i(T7|) are best possible in the 
sense that there exists a potential which turn these inequalities into equalities. This potential is such as W(r) = c, 
where c is a constant. The potential is then given by 

V(r)=a(f) -b(f) , (19) 

and c = AbR 2 . 



C. Upper limit for the lowest f-wave states 



To conclude this section devoted to the derivation of upper and lower limits on the rms radius, we present an upper 
limit applicable to the lowest £-wave states (no node in the wave function). This upper limit yields in general less 
restrictive constraints on the rms radius than the upper limits © and (|15p. except possibly for weakly bound systems. 

We begin with the inequality 



r 2 u 2 {r) 



< r 



dt(tu(t)Y 

OO 

dt [(tu(t))'} 2 



r / dtt 2 u'{ty, 



where the Cauchy-Schwarz inequality is used. The inequality Iplj) together with the relation J7J) yield 



-E(r 2 ) < 1+2 



drr 2 u'(r) 2 - 



where 



X 



1 + / drrV (r), 



(20) 



(21) 



(22) 



with V~(r) = max(0, — V(r)). The Jost-Pais necessary condition [24J implies that X is positive if the potential 
supports at least one bound state. The quantity X does not diverge if the potential decrease faster than r~ 2 at the 
infinity. Now, we consider the following inequality 



u'{r) 



dtu"(t), 
dt 



V(t) + V Z ' - E 



t 2 



< 



< 



dtV~(t)u(t) - 

i 

dtV-{t)u{t). 



dt 



E 



u(t), 

1(1+1) 
t 2 



u(t), 



(23) 
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The last inequality in Eq. ill'MI) is valid only for wave functions without node, i.e. the lowest £-wave states. The 
Cauchy-Schwarz inequality together with the inequality (|23|) yield 

/>oo />oc 

u' 2 (r) < / dt[V-{t)} 2 / dtu(t) 2 , 

J v J r 

< / dt[V-{t)] 2 . (24) 



We obviously suppose that the integral, from to oo, of the square of the negative part of the potential exists. This 
last inequality can be used with the relation (|21|l . and after an integration by parts we obtain 

-E(r 2 )<l + - drr 3 [V-{r)} 2 ~ £(£+!), (25) 
Jo 

where X is defined by 121'P . The right hand side of the inequality (|25f) behaves as third power of the strength of the 
potential. Clearly, for large values of the strength this upper limit yields poor constraint on the rms radius but for 
values of the strength close to the critical value (value at which a first £-wa.ve bound state appears), and if this critical 
value is small enough, we show in sec. |E|that this upper limit is better than those reported in sees. Inland [OBI 

III. CRITERIA FOR OCCURRENCE OF HALOS 

In this section, we use the lower limit l|18|) to obtain a criteria for the occurrence of S-wave quantum halo states. 
These states are threshold phenomena characterized b y la rge mean square radii and small binding energies. They 
occur in nuclear physics as halo nuclei (see for example |25l l26t I27L Y28L |29| ) and in molecular physics as weakly bound 



dimers (see for instance |30|, |3l|, |32( ) . For quantum halo states, the separation energy of the two bodies of the system 
(a nucleus and a nucleon or two atoms) is much smaller than the mean binding energy of the particles which compose 
eventually these bodies. These many-body systems can then be treated as two-body systems interacting through a 
potential. 

The idea of the criteria is the following: knowing the two bodies which compose the system as well as the central 
interaction, we determine for which binding energies quantum halo states exist. The application pres ented in this 
section is somewhat complementary to previous studies found in the literature (see for example fssl l34l Issl Isfj l37l| 1 . 

Quantum halo states are characterized by an extension far out into the classical forbidden region. We consider that 
a state with a binding energy E is a halo state if its rms radius is larger than its classical radius 

(r 2 ) 1/2 >*r , (26) 

with E = V(ro). The value of a can be estimated if we consider that in a quantum halo states the probability to 
find the particles with a interdistance greater than the classical interdistance is greater than 50% (see for example 
[27l l28l ] ) . The value of a obtained with this definition varies from 1.37 for the square well to 1.68 for a potential which 
decrease as r -3 at the infinity. In this work we take a = 2 in all numerical calculations for simplicity. 
With the lower limit (|18|l and the constraint (|26|l we can write 

(r 2 )> — = l —>a 2 rl. (27) 

V ' ~ -IE -2V(r ) " v ' 

Thus we are in the halo regime if — 2a 2 r^V (r^) < 1. 

We can now obtain the criteria for the occurrence of halo states. We write the potential in the convenient form 
V(r) — —gR~ 2 v(r / R). Notice that the number of bound states in the potential and the critical value of the strength 
g at which new bound states appear do not depend on R. We first search for the largest solution, xq, of the equation 

x 2 v(x ) = 1 , (28a) 
z<7 9c 

where xq = r /R and where the coupling constant g is replaced by its critical value g^ 1 for which the Nth eigenstate 
has a vanishing binding energy. Since g > g^ , we obtain a value of xq slightly underestimated which lead to a slightly 
lower value of the energy Eh above which quantum halo states appears. This effect will be attenuated (suppressed 
in practice) if we use upper limits on g^ . Accurate upper and lower limits on the value of g^ and on 3^ r>1 , involving 
only the potential, can be found in the literature EE E3, E3, EM 113, E3, 11M 11^, E3| . 
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The energy Eh is then given by 



E H = -g?R-*v(xo) = -^Aj- (28b) 

2(ax ) 2 



Consequently, states characterized by a binding energy E larger than Eh are characterized by a rms radius satisfying 
the inequality ij^Bt - 

Relevant information about occurrence of quantum halo sates is obtained with the criteria l|28[l only if the energies 
at which these states appear is obtained with a reasonable accuracy. In other words, the inequality (|26() should be 
verified and reasonably close to saturation. Two tests are performed below. 

Several remarks are in order. 

• Since the inequalities g\ < g 2 < . . . < g^ are always verified, x$ is larger for an excited state than for the 
ground state (see (JSSaJ). Consequently, the energy Eh is greater for excited states than for ground states. This 
clearly indicates that halo states are likely to be ground states instead of excited states. This conclusion may 
be incorrect for potentials which vanish identically beyond a given radius x* — r* / R. In this case, xq could 
stay constant for all values of N if the radius x* is small enough. This is the case for a square-well potential as 
discussed below. 



• From the relation J28bj, halos states have best chance to appear in potentials with a small range R. Indeed, 
when R varies, the quantities g^v(xo) or xq remain unchanged. This result is simple to understand: both 
(r 2 ) 1 / 2 and tq scale like R, consequently, their ratio is independent of R but the energy E scales like R~ 2 . 
Consequently, the energy for which the system is characterized by a given value of the ratio (r 2 ) 1 / 2 /ro scales 
like R- 2 . 

• When the two constituents of the halo are characterized by a finite size, like atoms or nuclei, halo states have 
best chance to appear for small sizes and for small reduced masses of these constituents (if a repulsion exist for 
small interdistance) . 



To illustrate the last affirmation, we consider the following interaction 



V{r) = gR- 



i?\ 2(n_1) (R 



r \ r 



(29) 



The repulsive part of the potential takes roughly into account the Pauli repulsion and R is then linked with the sizes 
of the particles interacting through this potential. The attractive part describe various kind of interactions depending 
on the value of n. The interaction of a charge and an induced dipole correspond to n — 4; both n = 6 and n = 7 
correspond to van der Waals forces, of London and Casimir-Polder type, respectively. For n — 6, we have a Lennard- 
jones (10,6) potential. The particular form of the interaction l|29[l is simply chosen to allow analytical calculations 
and is not intended to describe physical systems accurately but instead to get insights for the gross characteristics of 
some physical systems. 

The function W(r) (see l|14|l ) is positive for n > 4 and the lower limit (|18|l and then the criteria H28JI can be used. 
An estimation of the critical value g^ of the strength g is obtained with the formula 



9? < (N«) 2 
which is applicable for n > 4 39]. This leads to 



dx y/ v~(x) 



Uo 



(30) 



g? < (2N(n ~ 2)f. (31) 



The quantity xq is obtained easily 



where 7 = l/(2a 2 g^) < 0.0078 which justify the expansion around 7 = 0. The square root in the relation l|32|) 
imposes g^ > 1/2. This is always verified since the simple lower bound g^ > 2N(n — 2) can be obtained with the 
Bargmann- Schwinger inequality H [H) . The relations jSST^ . (|3TT) and lead to 

(2 n+A a 2n N 2 (n-2) A )^ 1 
E » = -- 2,R 2 1 ' (33 ) 
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where the reduced mass of the system, fi, has been written explicitely. The energy Eh decreases when n grows giving 
more chance to halo states to exist. For fixed value of n, Eh increases with jx and R and also with N. 

We can calculate Eh with the formula (|r?H|> for n — 6 and R — 2.640 A 41]. This model could be used to describe 
roughly the helium dimer 4 He2. We find that Eh = —0.82 /ieV, while the experimental binding energy is found to 
be around —0.095 fieV |42j. In this model, the helium dimer is a halo state. Actually the value of R used in this 
simple calculation was adjusted for a Lennard- Jones (12,6) potential |^3. If this last potential is used, instead of a 
Lennard-Jones (10,6) potential, we find a slightly modified value Eh = —0.89 fieV. 

We can also study how the energy Eh is sensitive to the asymptotic behavior of the potential. We simply choose 
the potential 

" -ITW (34) 

For n > 4, the formula l|30fl can be used to obtain an upper bound on the critical coupling constant. We have 

N 

■ 9c - r 2 (F ¥) r 2 (i + I) ' (35) 

where r(x) is the Euler gamma function. To find the quantity xq, we need to solve 

^ = ^-1-^, (36) 

7 7 

with 7 = l/(2cr 2 . g f ). This leads to 

x Q = (2a 2 g?)^. (37) 
The energy at which halo states appear is given by 

E H = -R- 2 (2v 2 )^(g?)^, (38) 
= -R- 2 s(N,n). (39) 

Before to draw some conclusions about the influence of the asymptotic behavior of the potential on the energy Eh, 
we mention that, as expected, the rms radius satisfy the inequality i|26|) if the binding energy is larger than Eh- If 
the binding energy is chosen to be equal to Eh for N = 1, an exact numerical calculation shows that the value of the 
ratio of the rms radius over the classical radius is, for example, 2.44 for n = 5, 2.47 for n — 10, 2.44 for n = 20, 2.41 
for n = 50 and 2.40 for n — 100. The same calculation for N — 2 leads the following values for the same ratio: 2.29 
for n = 5, 2.38 for n = 10, 2.42 for n = 20, 2.40 for n = 50 and 2.40 for n = 100. These results indicate that even if 
the function W(r) computed with the potential (|34|l is partially negative for n > 6, the criteria stays applicable for 
large values of n since the inequality Q26[) is always verified. 

Fig. ^indicates clearly that halos could exist more easily when n is large since in this case the energy at which they 
appear, Eh, is smaller (keeping R constant obviously). It is also clear that halos have more chance to exist as ground 
states than as excited states. However this last conclusion is only partially correct since we have linin^oo s(N,n) — 
l/(2er 2 ) for all value of N. In this limit the potential ffi^i reduces to a square- well potential. This property can be 
verified with an exact calculation, i.e. for a square-well potential, the energy at which halo states appears is almost 
the same for the ground state and for excited states. For example, with the potential V(r) = — Vq exp(— r/R)8(R — r), 
halo states have even marginally more chance to exist as excited state than as ground states. It is also clear on 
the Fig. n that the ratio of the energy Eh for the ground states over the energy Eh for the first excited states, 
Eh(N — 1)/Eh{N = 2), decreases as n grows. The values of this ratio are for example: 4 for n — 4, 2 for n = 6 and 
1.41 for n = 10. 

The same study could be performed with the potential V(r) = —gR~ 2 exp(—(r/R) n ). The same qualitative 
behaviors of Eh than those obtained with the potential (|34|) are observed. The ratio Eh(N = 1)/Eh(N = 2) also 
decreases as n grows but the values of this ratio are now smaller: 2.15 for n = 1, 1.35 for n = 2 and 1.12 for n = 5. 

We can also study how the energy Eh is sensitive to the details of the repulsive potential near the origin. We 
simply choose the potential 



V(r)=gR- 2 [ - - - . (40) 
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FIG. 1: s(N,n) as a function of n for several values of N. 
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For this potential, the function W(r) is positive for n > 6. For n > 5, the formula l|30|l can be used to obtain an 
upper bound on the critical coupling constant. We obtain 



g? < 4N 2 7t(n - 4) 2 r 2 ^ U 10 » ' 1 



To find the quantity Xq, we need to solve 



V 2n 



7^ = 1-4"" 



1 



1, 



with 7 = l/(2cr 2 5 ^). This leads to 

.T = (2cr g c ) 7 . 
The energy at which halo states appear is given by 

1 



E 



H 



R- 



(2a 2 ) 2 ^' 



-R~H{N, n). 



(41) 

(42) 
(43) 

(44) 
(45) 



Fig. El indicates that existence of halo states is not very sensitive to the details of the repulsive potential. For the 
ground state, the energy Eh increases by a factor 2 for n going from 5 to 25 whereas in the previous case, see Fig. ^ 
the increase was by a factor 10 in the same interval of n. 

Application to nuclear halos is possible. For this purpose we consider a Wood-Saxon potential as interaction between 
the nucleon and the nucleus 



V(r) 



V 



l + exp((r-i?)/a) : 



(46) 



with V = ga~ 2 , R = tqA 1 / 3 , r Q = 1.27 fm and a = 0.67 fm [43J. For this potential, the function W(r) is (slightly) 
partially negative whereas the criteria rigorously apply only to potential yielding a function W(r) everywhere positive. 
In Table [J we report a comparison between exact results and results obtained with the formulas l|28|) which prove 
that the criteria is applicable. 

The critical coupling constant of the potential (|46|l is simply estimated with the WKB formula 



jv 



(N- 1/4)tt 



a + 2 arcsinhl 



(47) 
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FIG. 2: t(N, n) as a function of n for several values of N. 



TABLE I: Comparison between the energy Eh given by l|49|l with the exact energy _E e x at which at which (r 2 ) 1//2 = 2ro for 
several values of the atomic number A. The ratio (r 2 ) 1 ^ 2 /ro at energy Eh is also given. The energies are given in MeV. 



A 


— -Bex 


-E H 


(r 2 } l/2 
TO 


1 


0.30 


0.26 


2.11 


5 


0.14 


0.11 


2.20 


10 


0.11 


0.087 


2.23 


15 


0.10 


0.075 


2.25 


20 


0.097 


0.067 


2.27 


25 


0.089 


0.062 


2.28 


•50 


0.069 


0.047 


2.30 


100 


0.053 


0.035 


2.32 



with a — R/a = 1.9 A 1 / 3 . The quantity Xq is obtained with the Lagrange inversion formula 

x Q = a (a 2 /7- 1) 1/q , 



(48) 



with 7 = l/(2a 2 gl). We restrict the calculation to the case N = 1 for which the formula l|48|l is quite accurate. 
Taking into account the reduced mass of the system, we finally obtain 



E H = -1.6 h(A) A- 2/3 MeV, 



where 



h(A) 



A+l 



A (o?h-\f la 



(49) 



(50) 



For a square- well potential, the energy En scales almost exactly as A" 2 / 3 , for a Wood-Saxon potential there is a 
correction. Indeed, taking the limit n — > oo of the relation l|38|l . and taking into account the reduced mass of the 
system, we obtain for a square-well Eh = —1.6 [(A + 1)/A] A~ 2 / 3 MeV (in agreement with a previous result |34|L 
The function h(A) is slightly increasing (h(l) = 0.16 and ft(225) = 0.55) and the modulus of the energy Eh decreases 
slower than A~ 2 I 3 . A numerical fit of the formula l)49|) leads to Eh — —0.22 yV 2 / 5 MeV. This last expression depends 
obviously on the value of a. If instead of a — 2 we choose cr = 1.5, the fit becomes E H = -0.25 A~ 1/2 MeV. 
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In Table [I] to test the criteria l|28[l applied to nuclear halos, we compare the value of the energy Eh given by the 
formula H49fl with the exact energy E eyi at which (r 2 ) 1 / 2 = ar , with a = 2. We also give the exact value of the ratio 
(r 2 ) 1 / 2 /tq at energy Eh- This ratio is almost constant, it increases slowly with the atomic number A. In any case, 
the inequality (|26[) is satisfied. 

IV. TESTS 

In this section we propose to test of the various upper and lower limits reported in sec.[n]with two simple potentials 



and 



V 2 (r)=- 9 R- ^- r/R) . (52) 



TABLE II: Comparison between the exact value of the quantity —E(r 2 ) for the ground state and various upper and lower limits 
reported in sec. [2] for the potential Vi(r) and for £ = 0. 



-E(r 2 ) Eq. © Eq. JT^J Eq. JEJ Eq. fift 



1.35 


0.51418 


1.7143 


1.0063 


0.5 


1.1549 


1.4 


0.55325 


1.7407 


1.025 


0.5 


1.1825 


1.5 


0.61348 


1.7937 


1.0625 


0.5 


1.2460 


1.75 


0.74037 


1.9259 


1.1563 


0.5 


1.4592 


2 


0.85531 


2.0582 


1.25 


0.5 


1.7623 


3 


1.26746 


2.5873 


1.625 


0.5 


4.1773 


4 


1.63531 


3.1164 


2 


0.5 


9.2480 


5 


1.97341 


3.6455 


2.375 


0.5 


17.949 



TABLE III: Same as Table El but for I = 1. 



g -E{r 2 ) Eq. © Eq. JTHJ Eq. ^ 



6.945 


0.06045 


2.6746 


1.771 


46.942 


6.95 


0.09321 


2.6773 


1.7729 


47.045 


7 


0.23641 


2.7037 


1.7917 


48.142 


7.5 


0.74323 


2.9683 


1.9792 


59.982 


8 


1.0652 


3.2328 


2.1667 


73.583 


9 


1.5923 


3.7619 


2.5417 


106.55 


10 


2.0529 


4.291 


2.9167 


148.03 



In Tables ITT1 and Hill we compare the exact value of the quantity —E(r 2 ), for the potential V\{r) and for £ = and 
1 respectively, with the various upper and lower limits reported in sec. [H] As expected the upper limit l|15H yields the 
strongest restrictions and is always better than the formula (JSJ) by construction. However for £ = 1, the constraints 
are poor for states characterized by a weak binding energy. As explain in sec. Ill CI the upper limit l|25[) yields its best 
restrictions for the coupling constant g close to the critical value g c (1.3326 < g c < 1.3403 |T3ll40| ). Useful restrictions 
are obtained only if this critical value is small enough. For £ = 1, we have g c 6.94 (6.9221 < g c < 6.9535 0IU) 
and the relation (|25|) gives poor restrictions. If g c is smaller than 1, the upper limit l|25|) could be better than i|15[) 
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for g close enough to g c ; this is the case for the potential V(r) — —gR~ 3 rexp(—r/R). For the potential Vi(r), the 
lower limit (|17fl yields non trivial results only for £ = 0. In this case, since the infimum of W(r) is equal to 0, the 
lower limit is a constant and is quite restrictive for values of g close to the critical value g c . In general, this lower 
limit yields strong restrictions when the binding energy of the system is small. This explain the good accuracy of the 
criteria obtained in sec. IIIII 

In Tables llVl and Ivl we compare the exact value of the quantity —E(r 2 ), for the potential T^(r) and for £ = and 
1 respectively, with the various upper and lower limits reported in sec. [n] For this potential, the upper limit (|25l) is 
not applicable and the lower limit l|17fl yields non trivial results only for £ = 0. Again, as expected the upper limit 
(I15|l yields the strongest restrictions and is obviously always better than the formula JSJ . The lower limit (|17fl is again 
very accurate for weak binding energy. 

TABLE IV: Same as Table EH but for the potential Va(r) and £ = 0. 



g -E{r 2 ) Eq. © Eq. Eq. JTJJ 



1.7 


0.50948 


1.6254 


0.92778 


0.49575 


1.75 


0.53314 


1.6438 


0.94036 


0.49563 


1.8 


0.55699 


1.6622 


0.95294 


0.4955 


1.9 


0.60473 


1.699 


0.97810 


0.49525 


2 


0.65199 


1.7358 


1.0033 


0.495 


3 


1.0558 


2.1036 


1.2549 


0.4925 


4 


1.3429 


2.4715 


1.5065 


0.49 


5 


1.55389 


2.8394 


1.7582 


0.4875 



TABLE V: Same as Table EH but for the potential V 2 (r) and I = 1. 



g -E{r 2 ) Eq. © Eq. 



9.085 


0.03585 


2.3422 


1.4528 


9.1 


0.092294 


2.3477 


1.4565 


9.5 


0.53225 


2.4949 


1.5572 


10 


0.85248 


2.6788 


1.683 


11 


1.3309 


3.0467 


1.9346 


12 


1.7085 


3.4146 


2.1863 


13 


2.0267 


3.7824 


2.4379 


14 


2.3031 


4.1503 


2.6895 


15 


2.5476 


4.5182 


2.9412 



V. CONCLUSIONS 



In this article, several rigorous upper and lower limits on the rms radius have been obtained for systems governed 
by central potentials. Some of these limits are applicable to eigenstates with arbitrary radial quantum number and 
angular momentum. Some of these limits yield in general strong restrictions on the rms radius as shown in sec. IIVI 
The simple lower limit l|18|) is used to obtain a criteria for the occurrence of S-wave halo states. This criteria gives the 
binding energy, called Eh , above which the eigenstate is characterized by a large rms radius compared to the classical 
radius, (r 2 ) 1 / 2 > ar$, and is thus qualified as quantum halo states. The relevance of the criteria has been tested with 
various potentials and we found that accurate information are obtained. It is worth noting that the various formula 
derived in sec. IIIII are applicable for arbitrary values of a but the various numerical values reported in the text and 
in the tables are computed for a — 2 for simplicity. This value could be changed for practical uses. However, the 
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conclusions obtained in this work, and summarized below, do not depend on the precise value of a contrary, of course, 
to the energy E H . 

With this criteria we have shown that halo states are likely to be ground states and not radial excitations (except 
possibly for potentials which vanish identically beyond a given radius, like the square-well potential, for which the 
converse could be true). This conclusion complete others results obtained previously which proved that halo states 
are S- or P-wave states. We have shown that, when the two constituents of the halo are characterized by a finite 
size, like atoms or nuclei, halo states have best chance to appear for small sizes and for small reduced masses of these 
constituents (if a repulsion exist for small interdistance). The criteria is also used to confirm that halo state have more 
chance to exist in potentials which tends rapidly to zero asymptotically. We have also shown that, if the potential 
has a repulsive part near the origin, the existence of halo is not very sensitive to the detail of this repulsive part. 
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